New effective interactions in RMF theory with non-hnear terms and 
density-dependent meson-nucleon couphng * 

Wenhui Long,^'^ Jie Meng,^' ^' ^ Nguyen Van Giai,^ and Shan-Gui Zhou^' ^' 

School of Physics, Peking University, 100871 Beijing, China 
^Institut de Physique Nucleaire, CNRS-IN2P3, Universite Paris-Sud, 91406 Orsay, France 
'"^Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing, China 
Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator, 730000 Lanzhou, China 
^ Max- Planck- Institut fur Kernphysik, 69029 Heidelberg, Germany 



CO ■ New parameter sets for the Lagrangian density in the relativistic mean field (RMF) theory, PKl 

' with nonlinear a- and oj-meson self-coupling, PKIR with nonlinear a-, u- and p-meson self-coupling 

' and PKDD with the density-dependent meson-nucleon coupling, are proposed. They are able to 

, provide an excellent description not only for the properties of nuclear matter but also for the nuclei 

^ ' in and far from the valley of beta-stability. For the first time in the parametrization of the RMF 

Q ' Lagrangian density, the center-of-mass correction is treated by a microscopic way, which is essential 

, to unify the description of nuclei from light to heavy regions with one effective interaction. 

Q ■ PACS numbers: 21.60.-n, 21.30.Fe, 21.60.Jz, 21.10.Dr 
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I. INTRODUCTION 



In the past decade, the development of unstable nuclear beams [1, 2] has extended our knowledge of nuclear 
physics from the stable nuclei and those nearby to the unstable nuclei far from the stability line. Intense research 
, in this area shows that there exist lots of unexpected phenomena: strange nuclear structure like neutron halo (skin) 
• and proton halo (skin) [3-10], soft excitation modes [11, 12], the enhancement of fusion cross sections induced by 
the extended matter distributions [13, 14] etc.. With further developments, many other new features will be found. 
It also becomes very important to find a reliable theory and improve the reliability for predicting the properties of 
even more exotic nuclei out to the proton and neutron drip lines. 
^ , Relativistic mean field (RMF) [15, 16] theory has received wide attention because of its successful description of 

many nuclear phenomena during the past years. With a very limited number of parameters, RMF theory is able 
to give a satisfactory description for the ground state properties of spherical [17] and deformed nuclei [18] at and 
away from the stability line. The recent reviews on RMF theory can be seen in [16-18]. In the simplest version 
of RMF theory, the mesons do not interact among themselves, which results in a too large incompressibility for 
nuclear matter. Boguta and Bodmer [19] therefore proposed to include a nonlinear self-coupling of the cr-field, 
' a concept which has been used in almost all the recent applications. The meson self-coupling introduces a new 
density dependence into the Lagrangian and consequently, the nuclear matter incompressibility can be lowered to 
reasonable values. As an implement, in 1994 the nonlinear self-coupling of the w-field is introduced by Sugahara 
and Toki [20]. In this paper we will introduce the nonlinear self-coupling for the p-field. Recently RMF theory 
with density-dependent (DD) meson-nucleon couplings [21-24] was developed by various authors. 

Till now the two versions (the nonlinear self-coupling of meson fields and the DD meson-nucleon couplings) of 
RMF theory have been successfully applied to describe the nuclear properties, including binding energies, nuclear 
matter distribution, single-particle spectra, magnetic moments, collective excited states, dipole sum rule, shell 
effects in finite nuclei, pseudospin symmetry, rotating nuclei, superdeformed bands, etc. In particular, the halo 
phenomena can be understood self-consistently in this microscopical model after the proper treatment of the 
continuum effect [7, 25]. Combining with Glauber model, the charge-changing cross sections for C, N, O and F 
isotopes are calculated and good agreement with the data has been achieved[9, 10]. The long existing problem for 
the origin of pseudo-spin symmetry in nuclei is given naturally as a relativistic symmetry [26-28] . Good agreement 
with experimental data has also been found recently for magnetic rotation [29] , collective excitations such as giant 
resonances [30] and for twin bands in rotating super-deformed nuclei [31]. It is also noted that cranked RMF 
theory provides an excellent description of superdeformed rotational bands in the A=140-150 region [32], in the Sr 
region [33] and in the Hg region [34]. 
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Among the existing parametrizations for RMF theory, the most frequently used are NLl [35], PL-40 [36], NL- 
SH [37], TMl [20] and NL3 [38] with nonhnear self-couphng of mesons, and TW99 [23] and DD-MEl [24] with DD 
mcson-nuclcon coupling. The effective interactions NLl, NL3, TM1,NL-SH and TW99, DD-MEl give good results 
in most of the cases. 

Along the beta stability line NLl gives excellent results for binding energies and charge radii, in addition it 

provides an excellent description of the superdeformed bands [32, 33]. However, in going away from the stability 
line the results are less satisfactory. This can be partly attributed to the large asymmetry energy J ~ 44 MeV 
predicted by this force. In addition, the calculated neutron skin thickness shows systematic deviations from the 
experimental values for the set NLl. 

In the effective interaction NL-SH this problem was treated in a better way and the improved isovector properties 
have been obtained with an asymmetry energy of J ~ 36 McV. Moreover, NL-SH seems to describe the deformation 
properties in a better way than NLl. However, the NL-SH parametrization produces a slight over-binding along 
the line of beta-stability and in addition it fails to reproduce successfully the superdeformed minima in Hg-isotopes 
in constrained calculations for the energy landscape. A remarkable disagreement between the two parametrizations 
are the quite different values predicted for the nuclear matter incompressibility. NLl predicts a small value {K 
= 212 MeV) while with NL-SH a very large value {K = 355 MeV) is obtained. Both forces fail to reproduce the 
experimental values for the isoscalar giant monopolc resonances for Pb and Zr nuclei. The NLl parametrization 
underestimates the empirical data by about 2 MeV while NL-SH overestimates it by about the similar values. 
As an improvement, the effective interactions, NL3 and TMl, provide reasonable compression modulus (i^NLs = 
268.0 MeV, Ktmi = 281.16 MeV) and asymmetry energy (Jnls = 36.56 MeV, Jtmi = 36.89 MeV) but fairly small 
baryonic saturation density (pnl3 = 0.145, pxMi = 0.145). In order to improve the description of these quantities, 
we developed two nonlinear self-coupling parametrizations called as PKl with nonlinear <t- and w-meson nonlinear 
self-coupling, and PKIR with nonlinear a-, oj- and /t>-meson nonlinear self-coupling. 

RMF theory with DD meson-nucleon couplings [21, 22] is an alternative approach to the description of nuclear 
matter and finite nuclei as compared to the model with the nonlinear self-interactions of mesons. There exist 
two representative parametrizations, TW99 [23] and DD-MEl [24] for the density-dependence of meson-nucleon 
coupling. They are able to describe quantitatively properties of nuclear matter and finite nuclei with similar quality 
as the parametrizations of the nonlinear self-coupling. As a comparison with the nonlinear self- coupling, we also 
developed a parametrization named as PKDD with DD meson-nucleon couplings. 

In all previous nonlinear self-coupling parametrizations as mentioned above, the center-of-mass correction is made 
by a phenomenological way. In the present parametrizations, the PK-series, the contribution from the center-of- 
mass motion is treated in a microscopic way [39]. The systematic behavior of the center-of-mass correction on 
nuclear masses is shown in Fig. 1. From this graph, we can see that it is essential to choose a proper method to 
treat the center-of-mass motion for both light and heavy nuclei. Obviously, the microscopic method provides more 
reasonable and reliable results for the center-of-mass motion (for the details, see Sec. III). 

In section II, we will present a short summary of RMF theory with the nonlinear self-coupling of meson fields 
and with DD meson-nucleon coupling and the relationship between them. The details of our parametrizations are 
given in section III. In section IV and V, we study the bulk properties of nuclear matter and spherical nuclei 
with the newly obtained effective interactions. The detailed microscopic structure of doubly magic nuclei is also 
investigated in the end of section V. Finally conclusions are given in section VI. 



II. RMF THEORY WITH NONLINEAR SELF-COUPLING AND DENSITY-DEPENDENT 

MESON-NUCLEON COUPLING 

The basic ansatz of the RMF theory [16] is a Lagrangian density whereby nucleons are described as Dirac 
particles which interact via the exchange of various mesons and the photon. The mesons considered are the scalar 
sigma (cr), vector omega (w) and iso- vector vector rho (p). The latter provides the necessary isospin asymmetry. 
The Lagrangian then consists of the free baryon and meson parts and the interaction part with minimal coupling, 
together with the nucleon mass M and rricr (ga), friui (Sw) and rUp {gp) the masses (coupling constants) of the 
respective mesons: 
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where the tensor quantities are 

a;'*'' = d^'io'' -d^iv^ (2a) 

p^'^ = df'p^ - + X (2b) 

Af" = d^" A" - d" Af" . (2c) 

In this paper we use arrow for isospin vectors and bold type for space vectors, respectively. There are 11 
parameters in the Lagrangian density (1), i.e., 4 masses {M,ma-,m^,mp), 3 nucleon-meson coupling constants 
{grr, Qlut gp), and 4 self-coupling constants (92, .93, C3, rfs). Generally, the nuclcon mass M and p-mcson mass m,p ( 
also sometimes the mass of w-meson in some parameterizations ) are fixed to their free values and the nonlinear 
coupling coefficients C3 and ds are taken as zeros. Then the remaining 6~7 parameters are determined by the 
fitting to the experimental observables. 

Different from the nonlinear self-coupling version, the meson-nucleon interactions arc described as density- 
dependent in the DD meson-nucleon coupling RMF theory [23, 24]. The meson-nucleon coupling constants g^, gw,gp 
become functions of the baryonic density py, p„ = ^/YpJ^ where = 'tp^y^tl; and the nonlinear self-coupling con- 
stants (?2j 53j C3, ^3 are set to zero in the Lagrangian density (1). 

For a and w-meson, the baryonic density-dependence of the coupling constants is adopted as 

gi {Pv ) = gi (Psat. ) fi {x) for i = a,uj (3) 

where 

l + hi{x + dif 

is a function of x = Pv/psat., and psat. denotes the baryonic saturation density of nuclear matter. 
For the p meson, an exponential dependence is utilized as 

gp = gp{Psa.t.)cxp[-ap{x - 1)] (5) 

For the functions fi{x), one has five constraint conditions /i(l) = 1,/^'(1) = /^'(l) and /-'(O) = 0. Then 
8 parameters related to density dependence for cr-N and w-N couplings are reduced to 3 free parameters. As 
mentioned above, the masses of nucleon and p-meson are fixed in general and the nonlinear self-coupling constants 
(/2,53,C3 and c?3 are set to zero. With 4 free parameters for density dependence, there totally arc 8^9 parameters 
left free in the Lagrangian density (1) for the density-dependent meson-nucleon coupling RMF theory. 

The single- nucleon Dirac equation is derived by the variation of the Lagrangian density (1) with respect to '0, 

[zj^dp ^ (M + Ss) - 7^S^] iP = (6) 
with the nucleon self-energies and E5 defined by the following relations: 

= ga(T (7a) 
^p = gu^'^p+gpT-pp + e^—^A^ + Y.^ (7b) 
where the rearrangement term comes from the density-dependence of the meson-nucleon coupling constants, 

= ^ {^i^l^^uj^ + ^ViT'-rV • p. + (8) 

Py \dpy dpy 8 Py J 

which is reduced to zero in RMF theory with the nonlinear self-coupling. 

The Klein-Gordon equations for mesons are obtained by the variation of the Lagrangian density (1) with respect 
to the corresponding meson field operators 

[-A -I- ma-] a = -gaPs - g2cr^ - gscr^ (9a) 
[-A -I- mj\ uj = gu,pb - csu^ (9b) 

[-A + mp]p = gp - - d^p' (9c) 

For the RMF theory with the nonlinear self-coupling, the self-coupling of mesons can be expressed into the density- 
dependence of meson-nucleon coupling by redefining the coupling constants, 

gcr = ga+ [ff2(T^ + 530-^] / Ps (10a) 

guj = gu>-c3uj^/pb (10b) 
gp = gp- dspV K - pI] (loc) 



4 



The behaviors of the couphng constants with respect to the baryonic density are shown in Fig. 2. 
In order to compare with the experimental mass of a nucleus, the calculated nuclear mass 

M = M„„f + Ec.ra. + £^pair (11) 

is obtained from the mean field contribution by adding corrections due to the center-of-mass motion and pairing 
effects for open shell nuclei. The mean field contribution Mmif is derived from a spatial integration of the "00" 
component of the energy-momentum tensor. In the nonlinear (N.L.) self-coupling case, it reads 



A 



-3„J 1 1 3 1 4 



r \ -^gcrcrps - -920- - -92.0 



(12) 



where (l>a denotes the single-particle spinor in the nucleus. In the density-dependent (D.D.) case, the high-order(> 3) 
terms of mesons in (12) are replaced by a "rearrangement" term, 

^rm?' =X^£a (^a7Va) + d^r ^^-^g^aps - ]^g^Ujpb - ^QpPpf^ " \eAoff^ - SflPfej . (13) 

where 

Sfl = — [5a.(Psat.)/<l,(a;)pow - ap5(p(psat.)e""''^'^"^Vo^V3 + 9a{ps^t.)f'a{x)psA (14) 

Psat. L J 

is the time-component of (see (8)). 

The correction from the center-of-mass motion is calculated from the projection-after-variation in first-order 
approximation [39]: 

E-^: = -^{K^) (15) 

where the center-of-mass momentum Pc.m. = X^f Pj and the expectation value of its square ^Pc.m.^ reads 

(Pc.m. ) = X] ~ S ^Wb'Pab • Po6 + X] ^a"«^b"i'Pa5 " Pa6 (16) 

a a,h a, 6 

with occupation probabilities and = 1 — accounting for pairing effects, where a, b denote the BCS states 
(see below). 

The contribution from the pairing correlations .Bpair are treated in the BCS approximation 

-Epair = - Ay^^tlgMg (17) 
a 

with the pairing gap A taken from the calculation of relativistic continuum Hartree-Bogoliubov (RCHB) theory 
with zero-range pairing interaction [25] . 

For the nuclear radii, the effects from the center-of-mass motion are also taken into account as follows. Because 
of its fairly small effects, a rather rough correction is adopted for protons 

A Z A 

SK = -^Y1 (-^"l • E + E (-^^l ^c.m. \^a) (18) 

a i a 

where the center-of-mass coordinate Rem. = iMX^f ^i- Then we get 

5Rl = -^Rl+^Rl, (19) 

where Rp and Rm denote the proton and matter radii. Here we only consider the direct-term contributions in 
(18) to keep with the spirit of RMF theory. For the neutron radii, we use the same procedure as for protons. The 
charge radius is obtained from the proton radius combining with the proton and neutron size and the center-of-mass 
correction (19) is included in R^ [20] 

^ch = Rl + (0-862 baf - (0.336 imfN/Z (20) 
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III. PARAMETRIZATION OF EFFECTIVE LAGRANGIAN AND NUMERICAL DETAILS 

The; aim of the present investigation is to provide new improved effective interactions for the Lagrangian density 
(1) with the nonUnear self-coupUng and density-dependent meson-nucleon coupling in RMF theory. A multi- 
parameter fitting was performed with the Levenberg-Marquardt method [40] . Two nonlinear self-coupling effective 
interactions have been obtained, PKl with a- and cj-mcson self-coupling, and PKIR with cr-, lo- and p-mcson self- 
coupling (see Table I). A density-dependent meson-nucleon coupling effective interaction has also been obtained 
with PKDD (see Table I and II). In our parametrization, the masses of neutron and proton are fixed to their free 
values: M„ = 939.5731 McV, Aip = 938.2796 McV, and the mass of p-meson is fixed to its experimental value 
763.0 MeV. The mass of w-meson is slightly adjusted in obtaining the effective interactions PKl and PKIR while 
fixed to 783.0 MeV for the effective interaction PKDD. 

The masses of the spherical nuclei ^^O., ""Ca, "^Ca, ^^Ni, ^SNi, ^"Zr, i^^Sn, 1323^^ i94pb 208pb are fitted to 
give the effective interactions PKl, PKIR and PKDD. The experimental inputs for finite nuclei used in the fitting 
procedure are shown in Table III. These inputs have been used to minimize the least square error: 



where a is the ensemble of parameters to be fitted, jy,^^^ and Gi are the experimental observable and corresponding 
weight. 

In the fitting procedure, only the masses of the spherical nuclei mentioned above and the compression modulus 
the baryonic density at saturation psat. and the asymmetry energy J of nuclear matter are included. Here, 
we should mention that the radii are excluded because we found that the values of the compression modulus and 
baryonic saturation density are essential to give a good description of the masses and radii. For a fixed value 
of compression modulus, a large baryonic saturation density will give a small charge radius. By this way, we 
can choose the appropriate value for these two quantities to give a proper description for both mass and charge 
radius of the chosen nuclei. To give a fairly precise description on the masses, the center-of-mass correction is 
essential for both light and heavy nuclei. As it can be seen in Fig. 1, the deviation between the microscopic and 
phenomenological results is considerably large not only for the light nuclei but also for the heavy ones. We also 
find that there exist very remarkable shell effects in the microscopic results which are impossible to obtain with 
the phenomenological methods. Hence, we choose the microscopic center-of-mass correction [39] to deal with the 
center-of-mass motion. 

The numerical procedure in obtaining the parameter set PKl is as follows. 

1. First, we start from an initial effective interaction and fix the mass of w-meson, let the other parameters to 
be adjusted by the Levenberg-Marquardt method [40]. 

2. Basing on the binding energies and charge radii of the selected nuclei obtained in the first step, we do 

some adjustment on the expectation of compression modulus and baryonic saturation density to improve the 
description on the spherical nuclei. Also, the weights of the observables would be slightly adjusted to improve 
the parametrization. 

3. Choose another initial effective interaction and take the same procedure as in step 2. Then, taking the average 
of these two obtained interactions as the new initial effective interaction, we do the minimizing procedure 
under the new weights which come from the previous results. 

4. Introduce the adjustment on the mass of cj-mcson to obtain the effective interaction PKl. 

Because the contribution to the nuclear mass from the nonlinear p-meson term is fairly small, we fix the nonlinear 
self-coupling constant rfs to 350.0 and adjust other parameters to obtain the effective interaction PKIR. 

To obtain the density-dependent meson-nucleon coupling effective interaction PKDD, we leave the density- 
dependent parameters a^, and d^^ to be adjusted and the five others to be determined by the constrain conditions 
on the functions fi{x),i = a,uj. The mass of w-meson is fixed to 783.0MeV. 

Thus obtained parameter sets are shown in Table I and II in comparison with other effective interactions TMl [20] , 
NL3 [38], TW99 [23] and DD-MEl [24]. In Table III and IV we respectively list the masses and charge radii of 
finite nuclei which have been used in obtaining these seven effective interactions. Compared with other effective 
interactions, our newly obtained ones reproduce well the experimental masses [41]. For these new effective inter- 
actions, only 4~5 nuclear massed deviate by more than IMeV (see Table III). The new effective interactions PKl, 
PKIR and PKDD also well describe the charge radii for these nuclei, especially for Pb isotopes. One can get a clear 
idea about the improvement of the new parameter sets on the description of bulk properties for finite nuclei from 
the root of relative square (rrs) deviation 5. In the last row of Table III (IV), the rrs deviation of the calculated 
total binding energy (charge radius) from the data is given. For the total binding energy, the rrs deviations from 
the new parameter sets are much smaller than those from old ones. For the charge radius, the rrs deviations from 




(21) 
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the new interactions are comparable with those from NL3 and DD-MEl, but a bit smaller than those from TMl 
and TW99. 

Tabic V lists the nuclear matter quantities calculated with the newly obtained effective interactions PKl, PKIR 
and PKDD, in comparison with other interactions. All the new effective interactions give proper value of the 
compression modulus K. Compared with the nonlinear self-coupling effective interactions TMl [20] and NL3 [38] , 
the new ones give more reasonable baryonic saturation density. 

IV. DESCRIPTION OF NUCLEAR MATTER 

We first discuss properties of nuclear matter obtained with PKl, PKIR and PKDD. In Table V we compare the 
bulk properties of nuclear matter with the results calculated using TMl, NL3, TW99 and DD-MEl. In Fig. 3, 
the behavior of the binding energy per particle E/A as a function of the baryonic density p is shown. We can 
see that all density-dependent mcson-nuclcon coupling effective interactions give softer results than the nonlinear 
self-coupling ones. The behaviors predicted by PKl, PKIR are much softer than by NL3 and a little harder than by 
TMl. The results from PKDD are slightly softer than that from DD-MEl and much harder than that from TW99 
at high density. All these behaviors can be explained in the density-dependent mcson-nuclcon coupling framework. 

As we have mentioned in expressions (10), the meson-nucleon coupling constants in the nonlinear self-coupling 
of mesons can be expressed as some kind of density-dependence. Fig. 2 shows this density-dependence for the 
nonlinear self-coupling effective interactions and the results for the density-dependent version are also given for 
comparison. We can see that almost all the density-dependent coupling constants decrease with increasing density 
except for g„ of NL3, which has a strong a self-coupling(g3 = —28.8851). On the other hand, the coupling constants 
Qcr and of TMl, which has relatively weak a self-coupling (33 = 0.6183) and strong uj self-coupling (cs = 71.3075), 
are smaller than the others, which means that TMl provides relatively weaker scalar and vector potentials. This 
is the reason why TMl presents the softer behavior than other nonlinear self-coupling effective interactions. In 
Fig. 3, TW99 predicts the softest results because of its relatively small as compared with DD-MEl, PKDD and 
NL3, large as compared with PKl, PKIR and TMl in Fig. 2. As we know, the repulsive potential would be 
dominant at high density. In Fig. 3, NL3 gives the hardest results because of its constant and large even though 
its a-N coupling constant g^ increase with the density. For the new parameter sets PKl, PKIR and PKDD, which 
present the mid soft behaviors in Fig. 3, the coupling constants also lie between the largest and the smallest in 
Fig. 2. For the parameter set PKIR, the density-dependence of the gp is fairly weak as compared with that of the 
density-dependent meson-nucleon coupling effective interactions. It can be explained by a very weak p-field, which 
generate neutron-proton symmetry field. The behavior of gp with respect to the neutron-proton ratio is shown in 
Fig. 4. As one can expect, the behavior is symmetric with respect to la.[N/Z) and the density-dependence becomes 
more remarkable with the increase of the baryonic density and the neutron-proton asymmetry. 

V. DESCRIPTION OF SPHERICAL NUCLEI 

A. Binding Energy and Two-neutron Separation Energy 

We calculate the even-even nuclei of Pb and Sn isotopic chains with the newly obtained parameter sets. In 
Table VI, we compare the masses calculated with PKl, PKIR, and PKDD with the other effective interactions and 
with experimental values [41]. Shown in Fig. 5 are the deviations of the masses of Pb isotopes from the data [41]. 
The results for Sn isotopes are shown in Fig. 6. In Fig. 5 and Fig. 6 we also give the results obtained with TMl, 
NL3, TW99 and DD-MEl for comparison. All results are calculated using the RCHB theory [25] where the pairing 
correlations are treated self-consistently by a zero-range ^-force. The box radius is 20 fm and the pairing strength 
of the zero-range (5-force is —650 MeV. The microscopic center-of-mass corrections (15) and (16) are used in all 
the calculations. 

As we can see in Table VI, Fig. 5 and Fig. 6, the newly obtained effective interactions PKl, PKIR and PKDD 
provide good descriptions on the masses of both two isotopic chains. In Fig. 5, all the effective interactions 
overestimate the binding energy in the beginning of the isotopic chain. However, from -"^^"Pb to ^^°Pb, the newly 
obtained interactions give better descriptions than all the others. For the Sn isotopes, the density-dependent 
effective interaction DD-MEl reproduces the data very well since 4 Sn isotopes were used in its parametrization [24]. 
From ^^^Sn to ^^^Sn, the new effective interactions PKl, PKIR and PKDD slightly underestimate the binding 
energy (less than 1 MeV). Compared with the others, the new ones still provide better description. There exist 
however systematic deviations out of the neutron magic numbers, e.g., in ^"^^Sn and ^^''Pb. For the older effective 
interactions, the deviations in ^^^Pb are smaller but fairly large in ^°^Pb. 

From the binding energies we can extract the systematics in two- neutron separation energies S2n = Eb {N, Z) — 
Eb{N — 2, Z). Figs. 7-9 exhibit two-neutron separation energies predicted by the nonlinear self-coupling effective 
interactions PKl, PKIR and the density-dependent meson-nucleon coupling one PKDD. In comparison, the results 
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obtained with TMl, NL3, TW99 and DD-MEl and the experimental values extracted from Ref. [41] are also given. 
In Fig. 10, the systematic behaviors of two-neutron separation energies with respect to neutron number, predicted 
by the new effective interactions PKl and PKDD, arc shown. From these figures one can sec that the newly 
obtained interactions give a fairly good description on the systematic behaviors in two-neutron separation energies. 
In Fig. 8, we can see that the deviations of theoretical results from experiment are rather large for Ni isotopes with 
N = 30-42 and for Sn isotopes with N = 52-58. Here it should be mentioned that all the theoretical results are 
extracted from the calculation of the spherical RCHB calculations [25] where deformation effects are not included 
while they can play a significant role in these nuclei. Furthermore, in Fig. 10 we can see some unusual phenomena 
along the neutron drip line which will be discussed in the next section. 

From the plots for two-neutron separation energies (Figs. 7-10), the position of the neutron-drip line for each 
element seems to be determined delicately. For the new effective interaction PKl, it predicts the neutron-drip 
number N = 50 for Ca, 70 for Ni, 96 for Zr and 126 for Sn. In general, the density-dependent meson-nucleon 
coupling effective interaction PKDD predicts smaller neutron number of the neutron-drip nucleus except for Ni. 
This may be due to its fairly small effective mass (see Table V) , which reduces the strong attractive potential in the 
core and makes the coupling between the core and valence orbital weaker. As for the density-dependent effective 
interactions, TW99 predicts smaller neutron numbers of neutron-drip nucleus for Ca, Ni and Zr while DD-MEl 
gives smaller ones for Ca and Sn. Another reason for the deviations among the effective interactions is that the 
density-dependent effective interactions give relatively large p-N coupling at lower densities (see Fig. 2). 

B. Charge Radius and Isotope Shift in Pb isotopes 

Although the radii are not included in our fitting procedure, the newly obtained effective interactions reproduce 
the charge radii of stable nuclei fairly well (see Table IV). The comparison between experimental data [42, 43] and 
theoretical results for the charge radii (20) of Pb isotopes are shown in Table VII. We can see that the new effective 
interactions PKl, PKIR and PKDD reproduce better the experimental values, as compared to overestimations by 
TMl, NL3 and DD-MEl and underestimations by TW99. We also calculate the isotope-shift of charge radii for 
Pb isotopes with these effective interactions. The results are shown in Fig. 11. The kink aromid ^"^Pb is well 
reproduced by all the interactions. The inset of Fig. 11 shows that the density-dependent meson-nucleon coupling 
effective interactions represent more reasonable agreement with the experimental values than the nonlinear self- 
coupling ones. The new parameter sets PKl and PKIR give slightly better results than TMl and NL3. 

Besides the charge radii, we have also used the new effective interactions to investigate the systematic behavior 
of the neutron skin (the difference between neutron and proton radii) along the isotope chains: O, Ca, Ni, Zr, Sn 
and Pb. Fig. 12 shows the radii difference r„ — rp calculated with the new sets PKl and PKDD. As we have seen in 
Fig. 10, there are several weakly-bound nuclei near the neutron drip line of Ca, Ni, Zr and Sn, whose two-neutron 
separation energies stay around zero over a range of numbers. This is in general a kind of signal for the existence 
of a neutron skin or halo. However from Fig. 12 one can see that the difference r„ — Vp tends to be a constant near 
the drip line for the Ni and Sn isotopes thus indicating a neutron skin rather than a neutron halo. On the other 
hand, the results support the existence of a neutron halo in Ca isotopes [44] and a giant neutron halo in Zr isotopes 
[8] because their neutron distributions tend to be more dispersive and rn — Vp keeps increasing rapidly. 

C. Single-Particle Energy and Spin-orbit Splitting 

The RMF theory is a microscopic theory with a limited number of parameters. It can give the detailed microscopic 
structure of nuclei. Figs. 13 to 16 show the single-particle energies of the doubly magic nuclei calculated with the 
newly obtained parameter sets PKl, PKIR and PKDD. In Figs. 13 and 14 the results obtained with TMl, NL3, 
TW99 and DD-MEl are also given for comparison. As we know, it is not straightforward to compare with the 
experimental results away from Fermi energies because of dynamical coupling not included in RMF. Here, the 
experimental values are extracted from one-nucleon separation energies [41] and resonance energies [45]. In these 
plots, we find systematic agreements with the experimental results. The single-particle energies near the magic 
numbers agree well with the experimental values. For levels far away from the Fermi energies systematic deviations 
appear. The states below the Fermi energy seems to be too strongly boimd, whereas the states above the Fermi 
energy show underbinding as compared to the experimental values. In fact, the experimental resonances are not 
simply single particle states. The coupling with core excitations through the residual interaction, which is obviously 
not included in the mean field model, leads in general to a shift of the resonance energies in the direction to the 
Fermi energy. In these plots of single-particle energies, the ordering of levels is well described by the new effective 
interactions except for somewhat different ordering of neutron levels in ^^^Sn and ^°^Pb, a common feature of all 
the effective interactions. 

From the single-particle energies mentioned above, we can extract the spin-orbit splittings. Table VIII shows 
the spin-orbit splitting calculated with the newly obtained effective interactions in the doubly magic nuclei. The 
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experimental values and the results calculated with other interactions are also given for comparison in Table VIII. 
The new effective interactions reproduce well the spin-orbit splittings. For the density-dependent meson-nucleon 
coupling effective interaction PKDD, the spin-orbit splitting turns out to be larger than the experimental values 
(except for the neutron 2p in '^^Ca and the neutron 3p in ^°*Pb). The nonlinear self-coupling effective interactions 
PKl and PKIR predict smaller splitting and improve the agreement with the experiment. The behavior can 
be explained by the relatively large self-energies predicted in our parametrization, which determine the strength 
of spin-orbit splitting. As shown in Table V, the nonlinear self-coupling effective interactions predict a larger 
effective mass than the density-dependent meson-nucleon coupling ones: the smaller the effective mass, the larger 
the spin-orbit splitting. This is also the reason why TMl gives the smallest splitting and TW99 the largest. 

VI. SUMMARY 

In this work, we have searched for new effective interactions to describe both stable and unstable nuclei in the 
relativistic mean field (RMF) theory with nonlinear self-coupling or density-dependent meson-nucleon coupling. 
In order to give a more precise description on the mass of nuclei, the microscopic center-of-mass correction is 
introduced, which makes it possible to give a unified description with one effective interaction for the nuclei from 
the light area to heavy area. As an elicitation from the density-dependent meson-nucleon coupling RMF theory, we 
introduce the nonlinear self-coupling for p-field. We obtain three new effective interactions: PKl with nonlinear 
self-coupling of cr-field and a;-field, PKIR with nonlinear self-coupling of cr-field, w-field and p-field and PKDD with 
density-dependent meson-nucleon coupling. 

With the newly obtained parameter sets, we investigate the behavior of the binding energy per particle and 
the meson-nucleon coupling constants with respect to the baryonic density in nuclear matter. The new sets PKl, 
PKIR and PKDD provide an appropriate description. Compared with TMl and NL3, the new ones give a more 
reasonable baryonic: saturation density. 

We then calculated the usual reference nuclei and Pb, Sn isotope chains and compared the masses with the 
available data. The nonlinear self-coupling effective interactions PKl, PKIR and density-dependent meson-nucleon 
coupling effective interaction PKDD reproduce well the data. As compared with other existing effective interactions, 
the new ones also provide a good description of the charge radii of the usual stable nuclei and Pb isotopes. 

We have also studied the systematics of two- neutron separation energies and neutron skin in isotopic chains. The 
two-neutron separation energies provided by the new interactions PKl, PKIR and PKDD are in good agreement 
with the experiment. We have also investigated the behavior of two-neutron separation energies and neutron skin 
near the neutron-drip line and given a reasonable interpretation of the formation of neutron halos. 

The single-particle energies and spin-orbit splittings in doubly magic nuclei predicted by the new parameter sets 
are compared with the experimental values and with other effective interactions. The new effective interactions PKl, 
PKIR and PKDD give a reasonable description of spin-orbit splittings and single-particle energies as compared 
with the experimental results. The systematic behavior of the spin-orbit splitting is interpreted in comparison with 
other effective interactions. 

Combining with the above information, we come to a conclusion that the new parameter sets PKl, PKIR and 
PKDD give better description for finite nuclei than other effective interactions. 
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TABLE I; The nonlinear effective interactions PKl, PKIR and density-dependent effective interactions PKDD. The masses 
(in MeV) and meson-nucleon couplings are shown in comparison with TMl [20], NL3 [38] and TW99 [23] and DD-MEl [24] 





PKl 


PKIR 


PKDD 


TMl 


NL3 


TW99 


DD-MEl 




939.5731 


939.5731 


939.5731 


938 


939 


939 


938.5000 


Mp 


938.2796 


938.2796 


938.2796 


938 


939 


939 


938.5000 


rria 


514.0891 


514.0873 


555.5112 


511.198 


508.1941 


550 


549.5255 


ruu, 


784.254 


784.2223 


783 


783 


782.501 


783 


783.0000 


rup 


763 


763 


763 


770 


763 


763 


763.0000 




1 n '^999 


in '^91 Q 




1 n n98Q 


1 n 91 RQ 


lu. / zoo 


1 n AA'\A 




13.0131 


13.0134 


13.1476 


12.6139 


12.8675 


13.2902 


12.8939 


9p 


4.5297 


4.55 


4.2998 


4.6322 


4.4744 


3.661 


3.8053 


92 


-8.1688 


-8.1562 





-7.2325 


-10.4307 





0.0000 


93 


-9.9976 


-10.1984 





0.6183 


-28.8851 





0.0000 


C3 


55.636 


54.4459 





71.3075 








0.0000 


ds 





350 

















TABLE II: Density-Dependent parameters of PKDD for meson-nucleon coupling in comparison with TW99 [23] and DD- 
MEl [24] 
















flp 


PKDD 


1.327423 


0.435126 


0.691666 0.694210 


1.342170 


0.371167 0.611397 


0.738376 


0.183305 


TW99 


1.365469 


0.226061 


0.409704 0.901995 


1.402488 


0.172577 0.344293 


0.983955 


0.515000 


DD-MEl 


1.3854 


0.9781 


1.5342 0.4661 


1.3879 


0.8525 1.3566 


0.4957 


0.5008 
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TABLE III: Total binding energies (in MeV) calculated with the nonlinear effective interactions PKl, PKIR and density- 
dependent meson- nucleon coupling effective interaction PKDD are shown in comparison with those of TMl [20], NL3 [38] 
, TW99 [23], DD-MEl [24] and experimental data [41]. The bold-faced quantities denote the observables used in the 
parametrizations. 



Nucleus 


Exp. 


PKl 


PKIR 


PKDD 


TMl 


NL3 


TW99 


DD-MEl 


16 Q 


-127.619 


128.094 


-128.047 


-127.808 


-128.951 


-127.127 


-128.147 


-127.926 


24q 


-168.500 


-169.558 


-169.381 


-168.542 


-168.858 


-170.116 


-167.693 


-167.949 


40Ca 


-342.052 


-342.773 


-342.741 


-342.579 


-344.661 


-341.709 


-343.352 


-343.653 


*«Ca 


-415.991 


-416.077 


-415.974 


-415.944 


-415.668 


-415.377 


-416.888 


-415.012 




-483.998 


-483.956 


-484.031 


-484.479 


-480.620 


-483.599 


-487.096 


-480.869 


58 Ni 


-506.454 


-504.033 


-504.091 


-504.013 


-501.933 


-503.395 


-506.128 


-501.312 


68 Ni 


-590.430 


-591.685 


-591.559 


-591.241 


-591.845 


-591.456 


-592.676 


-592.253 


90Zr 


-783.893 


-784.781 


-784.788 


-784.879 


-785.281 


-783.859 


-786.625 


-784.206 


"^Sn 


-953.529 


-954.210 


-954.251 


-953.730 


-955.925 


-952.562 


-954.991 


-952.468 


"«Sn 


-988.681 


-988.491 


-988.460 


-988.066 


-990.083 


-987.699 


-989.842 


-988.470 


^^*Sn 


-1049.963 


-1049.162 


-1048.948 


-1048.113 


1049.832 


1049.884 


-1051.033 


1049.880 


i32Sn 


-1102.920 


1103.503 


1103.053 


-1102.648 


-1102.163 


1105.459 


-1108.363 


1103.857 


184p^ 


-1431.960 


-1435.548 


-1435.706 


-1435.477 


-1439.768 


-1434.569 


-1436.761 


-1434.569 


194p^, 


-1525.930 


-1525.536 


-1525.494 


-1525.474 


-1528.378 


-1525.733 


-1529.309 


-1524.937 


196pb 


-1543.250 


-1542.592 


-1542.502 


-1542.545 


-1545.163 


-1543.085 


-1547.012 


-1542.262 


204p^ 


-1607.520 


-1607.851 


-1607.545 


-1607.770 


-1609.477 


-1609.906 


-1608.246 


1609.676 


208pb 


-1636.446 


-1637.443 


-1637.024 


-1637.387 


-1638.777 


-1640.584 


1644.790 


1641.415 


214p^ 


-1663.298 


-1659.382 


-1658.718 


-1656.084 


-1663.706 


-1662.551 


-1656.086 


-1662.011 


210p^ 


-1645.228 


-1648.443 


-1648.102 


-1648.039 


-1650.819 


-1650.755 


-1654.271 


-1651.482 


A" 




7.1980 


7.4207 


9.2744 


12.8135 


9.1140 


17.6762 


11.1580 


5* 




0.0102 


0.0094 


0.0080 


0.0192 


0.0135 


0.0159 


0.0152 



"The total square deviation from the experimental values A'^ = (B^"''' — ^'f*' ) 
''The relative square deviation 5^ = Ei i^T'^' - E9^^-f /E^'^'^-^ 
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TABLE IV: Charge radii (in frn) calculated with the nonlinear effective interactions PKl, PKIR and density-dependent 
meson- nucleon coupling effective interaction PKDD are shown in comparison with those of TMl [20], NL3 [38] , TW99 [23], 
DD-MEl [24] and experimental data [42]. 



Nucleus 


Exp. 


PKl 


PKIR 


PKDD 


TMl 


NL3 


TW99 


DD-MEl 


160 


2.693 


2.6957 


2.6959 


2.6988 


2.7026 


2.7251 


2.6799 


2.7268 


24Q 




2.8106 


2.8108 


2.8184 


2.8364 


2.8286 


2.8049 


2.8543 


40Ca 


3.478 


3.4433 


3.4435 


3.4418 


3.4541 


3.4679 


3.4151 


3.4622 


«Ca 


3.479 


3.4675 


3.4675 


3.4716 


3.4911 


3.4846 


3.4510 


3.4946 






3.7085 


3.7084 


3.7162 


3.7471 


3.7122 


3.6867 


3.7315 




3.776 


3.7383 


3.7381 


3.7442 


3.7755 


3.7435 


3.7158 


3.7613 


68Nj 




3.8621 


3.8620 


3.8681 


3.8901 


3.8773 


3.8491 


3.8926 




4.270 


4.2522 


4.2521 


4.2534 


4.2799 


4.2689 


4.2278 


4.2725 


"^Sn 


4.593 


4.5704 


4.5701 


4.5722 


4.6021 


4.5861 


4.5461 


4.5901 


"'^Sn 


4.625 


4.5984 


4.5981 


4.6004 


4.6303 


4.6149 


4.5758 


4.6212 


i^^Sn 


4.677 


4.6536 


4.6533 


4.6567 


4.6874 


4.6685 


4.6331 


4.6781 


i^^Sn 




4.7064 


4.7061 


4.7102 


4.7442 


4.7183 


4.6842 


4.7270 


184p^ 




5.3806 


5.3801 


5.3807 


5.4156 


5.3996 


5.3511 


5.4002 


194p^ 


5.442 


5.4327 


5.4322 


5.4329 


5.4712 


5.4506 


5.4017 


5.4539 




5.449 


5.4438 


5.4433 


5.4440 


5.4826 


5.4614 


5.4123 


5.4645 


204pb 


5.482 


5.4869 


5.4864 


5.4877 


5.5261 


5.5027 


5.4837 


5.5038 


208pb 


5.504 


5.5048 


5.5043 


5.5053 


5.5444 


5.5204 


5.4750 


5.5224 


2"Pb 


5.559 


5.5658 


5.5653 


5.5635 


5.6052 


5.5820 


5.5603 


5.5779 


210p^ 




5.5370 


5.5365 


5.5371 


5.5762 


5.5539 


5.5070 


5.5544 


A" 




0.0708 


0.0712 


0.0655 


0.0941 


0.0625 


0.1436 


0.0588 


5' 




0.0178 


0.0179 


0.0166 


0.0185 


0.0169 


0.0346 


0.0163 



"The total square deviation from the experimental values = (''™^ ~ ,.cai.j 
''The relative square deviation 5^ = ^ . (rj''"- - r?^' )^ /rf"^'^ 



TABLE V: Nuclear matter properties calculated with the nonlinear effective interactions PKl, PKIR and the density- 
dependent effective interactions PKDD are shown in comparison with TMl [20], NL3 [38], TW99 [23], DD-MEl [24]. 



Interaction 


Psat.(fm ^) 


Eb [MeV] 


K [MeV] 


J [MeV] 


M'/Min) 


M7M(p) 


PKl 


0.148195 


-16.268 


282.644 


37.641 


0.605525 


0.604981 


PKIR 


0.148196 


-16.274 


283.674 


37.831 


0.605164 


0.604620 


PKDD 


0.149552 


-16.267 


262.181 


36.790 


0.571156 


0.570565 


NL3 


0.145115 


-16.005 


267.998 


36.558 


0.603761 


0.603761 


TMl 


0.145218 


-16.263 


281.161 


36.892 


0.634395 


0.634395 


TW99 


0.153004 


-16.247 


240.276 


32.767 


0.554913 


0.554913 


DD-MEl 


0.151962 


-16.201 


244.719 


33.065 


0.577960 


0.577960 
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TABLE VI: Total binding energies of Pb isotopes (in MeV) calculated with the nonlinear self-coupling effective interac- 
tions PKl, PKIR and the density-dependent meson-nucleon coupling effective interaction PKDD, in comparison with the 
experimental values [41] and the results of TMl, NL3, TW99 and DD-MEl. 



A 


Exp. 


PKl 


PKIR 


PKDD 


TMl 


NL3 


rxTXT J C\C\ 

1 W99 


DD-MEl 


182 


—1411.650 


—1416.431 


—1416.619 


— 1416.309 


— 1420.872 


— 1415.184 


— 1417.202 


— 1415.216 


184 


—1431.960 


—1435.548 


—1435.706 


—1435.477 


— 1439.768 


■1 A^ A /"* r\ 

—1434.569 


—1436.761 


■1 A ^ A ^ /"* r\ 

—1434.569 


186 


—1451.700 


—1454.258 


—1454.382 


— 1454.192 


— 1458.222 


— 1453.511 


-1455.879 


-1453.306 


188 


—1470.900 


—1472.586 


—1472.673 


— 1472.504 


— 1476.272 


-1472.056 


— 1474.640 


-1471.639 


190 


—1489.720 


—1490.555 


—1490.603 


— 1490.468 


— 1493.958 


— 1490.247 


-1493.109 


-1489.657 


192 


—1508.120 


—1508.197 


—1508.201 


—1508.116 


— 1511.317 


— 1508.130 


— 1511.329 


— 1507.411 


194 


—1525.930 


—1525.536 


—1525.494 


— 1525.474 


— 1528.378 


— 1525.733 


— 1529.309 


— 1524.937 


196 


—1543.250 


— 1542.592 


— 1542.502 


— 1542.545 


-1545.163 


-1543.085 


— 1547.012 


— 1542.262 


198 


— 1560.070 


— 1559.378 


— 1559.236 


-1559.329 


-1561.685 


-1560.199 


-1564.338 


-1559.403 


200 


-1576.365 


-1575.893 


-1575.697 


-1575.831 


-1577.942 


-1577.076 


-1581.344 


-1576.365 


202 


-1592.202 


-1592.095 


-1591.843 


-1592.022 


-1593.905 


-1593.679 


-1598.084 


-1593.133 


204 


-1607.520 


-1607.851 


-1607.545 


-1607.770 


-1609.477 


-1609.906 


-1614.197 


-1609.676 


206 


-1622.340 


-1623.126 


-1622.765 


-1623.167 


-1624.530 


-1625.725 


-1630.122 


-1625.966 


208 


-1636.446 


-1637.443 


-1637.024 


-1637.387 


-1638.777 


-1640.584 


-1644.790 


-1641.415 


210 


-1645.568 


-1644.844 


-1644.345 


-1643.643 


-1647.245 


-1647.969 


-1650.888 


-1648.312 


212 


-1654.524 


-1652.155 


-1651.573 


-1649.873 


-1655.549 


-1655.289 


-1657.020 


-1655.174 


214 


-1663.298 


-1659.382 


-1658.718 


-1656.084 


-1663.706 


-1662.551 


-1663.196 


-1662.011 



TABLE VII: Chaxge radii of Pb isotopes (in fm), calculated with the nonlinear effective interactions PKl, PKIR and the 
density-dependent meson-nucleon coupling effective interaction PKDD, in comparison with those of TMl, NL3, TW99, 
DD-MEl and experimental values [42, 43]. 



A 


Exp. 


PKl 


PKIR 


PKDD 


TMl 


NL3 


TW99 


DD-MEl 


190 


5.4273 


5.4112 


5.4107 


5.4115 


5.4486 


5.4295 


5.3814 


5.4325 


192 


5.4347 


5.4219 


5.4214 


5.4221 


5.4599 


5.4400 


5.3915 


5.4433 


194 


5.4416 


5.4327 


5.4322 


5.4329 


5.4712 


5.4506 


5.4017 


5.4539 


196 


5.4487 


5.4438 


5.4433 


5.4440 


5.4826 


5.4614 


5.4123 


5.4645 


198 


5.4564 


5.4551 


5.4546 


5.4555 


5.4940 


5.4722 


5.4241 


5.4748 


200 


5.4649 


5.4664 


5.4659 


5.4673 


5.5052 


5.4830 


5.4366 


5.4850 


202 


5.4741 


5.4773 


5.4768 


5.4786 


5.5161 


5.4934 


5.4493 


5.4947 


204 


5.4820 


5.4869 


5.4864 


5.4877 


5.5261 


5.5027 


5.4571 


5.5038 


206 


5.4930 


5.4954 


5.4949 


5.4960 


5.5352 


5.5109 


5.4653 


5.5125 


208 


5.5040 


5.5048 


5.5043 


5.5053 


5.5444 


5.5204 


5.4750 


5.5224 


210 


5.5231 


5.5252 


5.5247 


5.5247 


5.5645 


5.5411 


5.4935 


5.5409 


212 


5.5415 


5.5455 


5.5450 


5.5441 


5.5846 


5.5616 


5.5122 


5.5594 


214 


5.5591 


5.5658 


5.5653 


5.5635 


5.6052 


5.5820 


5.5310 


5.5779 
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TABLE VIII: The theoretical (calculated with the nonlinear effective interactions PKl, PKIR, TMl [20], NL3 [38] and the 
density-dependent ones PKDD, TW99 [23], DD-MEl [24]) and experimental spin-orbit splittings (in MeV) of neutron (i/) 
and proton (tt) levels in doubly magic nuclei. 



Nucleus 


State 


PKl 


PKIR 


PKDD 


Exp. 


TMl 


NL3 


TW99 


DD-MEl 


16o 


flp 


6.550 


6.550 


6.950 


6.180 


5.660 


6.480 


7.480 


6.320 




Trip 


6.490 


6.500 


6.900 


6.320 


5.610 


6.400 


7.410 


6.250 




uli 


7 All 


7.417 


8.028 


8.380 


6.486 


7.493 


8.687 


7.498 




u2p 


1.238 


1.237 


1.459 


2.020 


1.142 


1.330 


1.567 


1.458 




flf 


8.223 


8.231 


8.670 


7.160 


6.907 


8.703 


9.274 


8.067 
























i/2p 


1.141 


1.142 


1.442 


1.110 


1.105 


1.112 


1.573 


1.388 




v2d 


1.659 


1.662 


1.990 


1.650 


1.515 


1.661 


2.257 


1.940 




Trig 


5.900 


5.910 


6.450 


6.080 


5.010 


6.152 


7.108 


6.210 




7r2d 


1.704 


1.706 


2.005 


1.750 


1.556 


1.690 


2.203 


1.893 




i/2f 


2.005 


2.008 


2.356 


1.770 


1.812 


2.011 


2.648 


2.268 




z/li 


6.492 


6.503 


7.126 


5.840 


5.634 


6.665 


7.761 


6.748 




i>3p 


0.742 


0.742 


0.879 


0.900 


0.657 


0.764 


0.992 


0.866 




7r2d 


1.626 


1.622 


1.832 


1.330 


1.436 


1.628 


2.031 


1.736 




Trlh 


5.448 


5.458 


5.976 


5.560 


4.653 


5.661 


6.576 


5.749 
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FIG. 1: The microscopic center-of-mass correction, in comparison with two phenomenological cases. 




FIG. 2: The meson-nucleon coupling constants as a function of baryonic density in nuclear matter. The density-dependence 
of Qp for PKIR corresponds to non-symmetric nuclear matter (N/Z = 3). 
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FIG. 3: The binding energy per particle E j A in nuclear matter as function of the baryonic density p, calculated with the 
density-dependent meson-nucleon coupling effective interactions PKDD, TW99, DD-MEl and the nonlinear effective inter- 
actions PKl, PKIR, TMl, NL3. The shaded area indicates the empirical value and the filled circles represent corresponding 
saturation points. 




FIG. 4: The density-dependence of for PKIR with respect to the neutron-proton ratio N/Z. 
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Mass Number 

FIG. 5: The deviation of the theoretical binding energies of Pb isotopes, calculated with the nonlinear effective interactions 
PKI, PKIR and density-dependent rneson-nucleon coupling effective interactions PKDD from the experimental values[41]. 
The results calculated in TMl, NL3 and TW99, DD-MEl are shown for comparison. 
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FIG. 6: Same as Fig. 5, for Sn isotopes. 
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FIG. 7: The two-ncutron separation energies in O and Ca isotopes calculated with the nonlinear effective interactions PKl, 
PKIR and the density-dependent rneson-nucleon effective interaction PKDD, in comparison with those of TMl, NL3, TW99 
and DD-MEl and the experimental data[41]. 
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FIG. 8: Same as Fig. 7, for Ni and Zr isotopes. 
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FIG. 9: Same as Fig. 7, for Sn and Pb isotopes. 
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FIG. 10: The two-neutron separation energies calculated with the nonlinear effective interaction PKl and the density- 
dependent meson-nucleon coupling one PKDD, as a function of the neutron number. 
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FIG. 11: The isotope shifts of the charge radius in Pb isotopes (in fm^), calculated with the nonlinear effective interactions 
PKl, PKIR and the density-dependent nieson-nucloon coupling effective interaction PKDD, in comparison with those of 
TMl, NL3, TW99, DD-MEl and experimental values, where Ar? = rl{A) - r?(208) and Ari^ = rlo{A) - rlo{208). 




FIG. 12; The radius difference r„ — Vp calculated with the nonlinear effective interaction PKl and density-dependent 
meson-nucleon coupling one PKDD, with respect to the neutron number. 
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FIG. 13: The single-particle energies in Pb, calculated with PKl, PKIR and PKDD, in comparison with the results of 
TMl, NL3, TW99, DD-MEl and experimental values. 




FIG. 14: Same as Fig. 13, for i^zgn. 
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FIG. 15: The single-particle energies in ■^^Ca and '^^Ni, calculated with PKl, PKIR and PKDD, in comparison with 
experimental values. 
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FIG. 16: Same as Fig. 15, for i«0 and *°Ca. 



